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Abstract. The recently published spectroscopic metallicities of RR Lyrae stars in M3 ISandstrom et all200ll) though show a 
relatively wide range of the [Fe/H] values, the conclusion that no metallicity spread is real has been drawn, as no dependence on 
either minimum temperature or period was detected. Comparing these spectroscopic metallicities with [Fe/H] calculated from 
the Fourier parameters of the light curves of the variables a correlation between the [Fe/H] values appears. As a consequence of 
the independence of the spectroscopic and photometric metallicities, this correlation points to the reality of a metallicity spread. 
The absolute magnitudes of these stars follow a similar trend along both the spectroscopic and photometric metallicities as the 
general M v - [Fe/H] relation predicts, which strengthens that the detected metallicity range is real. 
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1. Introduction 

Nearly every globular cluster contains stars of homogeneous 
metallicity. The most prominent exception is a> Cen with its 
considerably large, ~1 dex spread in the [Fe/H] content of 
giants and subgiants. Slight [Fe/H] inhomogeneity in M92 
king et alJfl99l iLaneer eUul ll998) has also been detected. 
Concerning the other elemental abundances, detailed chemical 
composition analysis of globular cluster stars have revealed that 
complicate patterns of abundance inhomogeneity exist in many 
clusters. 

M3 is one of the reference clusters of the globular clus- 
• ter metallicity scale. Spectroscopically, [Fe/H] has been ai- 
re ady derived from high dispersion studies of about 10 giant 
stars JBell & Dickens! 1980llKraft et alJl992HCavallo & Nagad 
2000; Kraft & Ivans 2002) which do not show any evidence 
of a [Fe/H] spread. The latest result gives an average [Fe/H] 
value of —1.5 dex on a metallicity scale based on measure- 
ments of Fell lines in order to avoid the non-n egligible bias 
of overionization effects in Fel lines JKraft & Ivansl20 02). The 
first spectroscopic data on direct [Fe/H] measurements of the 
cluster RR Lyrae vari ables have been recently published by 
ISandstrom et alJfcOOlL hereafter SPS). Using moderate resolu- 
tion spectra taken with the Hydra multifiber spectrograph they 
conclude that the RR Lyrae variables are also uniform in com- 
position as no dependence on either minimum temperature or 
period can be found based on data of 29 RR Lyrae variables. 
Although the scatter in SPS's [Fe/H] values ([Fe/H] mea n = 
-1.43, s.dev. = +0.12 and [Fe/H] mean = -1.21, s.dev. = ±0.22 
from Fel and Fell lines, respectively) can be explained by un- 
certainties of model parameters and equivalent widths, it does 



not exclude the possibility of some real star to star [Fe/H] abun- 
dance diff erences either. 

In ljurcsik & Kovacsl |l996) it has been shown that the 
[Fe/H] of RRab stars calculated from the period (P) and the 
03i Fourier phase-difference of the light curve is as accurate as 
most of the spectroscopic results. With this formula the spec- 
troscopic [Fe/H] of the 8 1 calibrating RRab stars could be re- 
produced with 0.14 dex standard deviation. 

In a preliminary investigation it has been found that the 
metallicities of the RRab stars in M3 calculated according to 
this photometric method show larger spread than expected if 
the [F e/Hl of the variables is homogeneous feakos & Jurcsikl 
2000). 

The aim of the present Letter is to check the reality of any 
possible [Fe/H] inhomogeneity of the M3 variables comparing 
their spectroscopic and photometric metallicity values. 

2. Spectroscopic and photometric [Fe/H] 

The two metallicity determination methods - namely, a) spec- 
troscopic; b) utilizing light curve parameters (photometric) -, 
are completely independent, therefore any correlation between 
these metallicities would indicate that the observed metallicity 
spread is real. 

2.1. Spectroscopic metallicities 

During the last few years strong arguments were raised against 
metallicity determinations using Fel lines dLambert et al.ll 996 ; 
iThevenin & Idiartlll999t iKraft & Ivansll2002l) because of the 
non-negligible departures from LTE in metal-poor stellar at- 
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mospheres resulting overionization of iron. Based on these re- 
sults for the comparison purposes we use exclusively the Fell 
metallicities of SPS, which have unexpectedly large, > 0.6 dex 
range. 

The signal-to-noise ratio (S/N) of the spectroscopic obser- 
vations varies between 28 and 76 whereas the number of spec- 
tra (NS) of individual stars are between 10 and 41 . For V9 both 
the S/N and NS have one of the lowest values, 30 and 10, re- 
spectively. Therefore we have decided to omit this star from 
the spectroscopic sample. The [Fe/H]Feii = -0.87 dex seems 
to be unreliably large for this star based on the comparison of 
its light c urve with a la rge sample of RRab light curves com- 
piled in IJurcsikl dl998h . A field RR Lyrae star (UW Gru) has 
been found to show light changes very similar to M3/V9. 



value cannot be correct for M3/V9. It validates our decision to 
omit V9 from the spectroscopic sample. 

Table 1. Parameters of M3/V9 and UW Gru 



UWGru 



ref. 



M3/V9 ref. 



period [d] 
(B - V) [mag] 
E(B-V) [mag] 
[Fe/H] phot [dex] 



0.5482104 
0.36 
0.10 
-1.51 



[Fe/H] spec , [dex] -1.68,-1.41 5,6 



0.5415552 
0.33, 0.39 
0.01 
-1.53 
-0.87 



2.3 
4 



1) Bernard ( 1982); 2) Benko et al. 
|200l| ) ; 4) |HarrisHl99fj); 5 ) |Lavden 
1 1996); 7) Sandst romet all 1200 ll) 



Corwin & Carnev 



Jurcsik & Kovacs 




Fig. 1. Phased V light curves of M3/V9 (circles) and UW Gru 
(crosses). The magnitudes of UW Gru are shifted to match the 
magnitudes of M3/V9. The identity of the light curves' shapes 
indicates close similarity of the objects that have similar peri- 
ods and colours. Consequently, no significant metallicity dif- 
ference between these stars can be real. 



Fig. ^ shows the phased li ght curves of M3/V9 and UW 
Gru using the obs ervations of ICorwin & Carnevl J200ll) and 
iBenko et alJ d2003l) for M3/V9. and lBernardl Jl982l) for UW 
Gru. The light curves of these two stars can be hardly sepa- 
rated, within photometric uncertainties they can be regarded 
as identical. We suppose that such equality of the light curves 
with similar pulsation periods can only happen if there is no 
significant difference in any of the basic physical properties 
of the two stars. The periods of these variables are very sim- 
ilar and their mean colours are also the same within the range 
limited by the photometric uncertainties as given in Table [T] 
The metallicity of U W Gru was given as -1.68 dex in ( Lavden 



1994) which was transformed to -1.41 dex metallicity value 
on a common [Fe/H] scale used in Jurcsi k & Kovacsl {l996). 
As a consequence of the absolute equality of the light curves 
of M3/V9 and UW Gru the spectroscopic [Fe/H] = -0.87 dex 



2.2. Photometric metallicities 

Reliable photometric [Fe/H] can be determined for 19 RRab 
stars among the spectroscopic sample of SPS (29 stars). The 
following stars are omitted: 

- stars showing strong Blazhko modulation that makes 
[Fe/H] phot incorrect (V20, V39, V52, V61, V63, V66, V78). 

- no or poor quality published light curves are available 
(V113,V115). 

- V202, which cannot be regarded as a normal RRab star 
based on its long period ((fill), light curve shape, and low 
amplitude (A v = (T 18). 

In order to increase the small sample size those Blazhko stars 
are also used in the comparison which exhibit only slight light 
curve modulations (V10, V59 and V62). The small amplitude 
modulations are supposed not to affect significantly the calcu- 
lated photometric metallicities. 

For most of the stars photometries are available from dif- 
ferent sources. In these cases the data have been merged with 
necessary zero-point offsets (0.01 - 0.05 mag) and a single 
[Fe/H]ph ot is calculated from the combined data set. In each 
case the magnitud e correction has been determined relative 
to the ICorwin & Ca rney (2001) data. All the data have been 
cleaned from outlying points. The large number of the photo- 
metric measurements (N« 120 - 700) makes the separation of 
discrepant data unambiguous in most of the cases. The formal 
errors of the [Fe/H] p h t values are between 0.03 and 0.07 dex. 

2.3. [Fe/H] 

spect 

[Fe/H] ph 

ot 

The spectroscopic and photometric metallicit ies, the inten- 
sity a veraged mean magnitudes according to lCorwin & Carnevl 
(2001), and the used photometries of the variables are listed in 
Tabled 

The final sample of stars with [Fe/H] reliably determined 
both spectroscopically and photometrically consists of a lim- 
ited number of 18 variables. Fig. |2] shows these data, with 
the direct and inverse regression fits overplotted. The linear 
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Table 2. Spectroscopic and photometric metallicities and ob- 
served mean magnitudes of the SPS RRab sample. 



Var. 


[Fe/H]* 


[Fe/H] phot 




Rem. b 


Ref. phot. c 




dex 


mag 






9 


-0.87 


-1.53 


15.630 


* 


2, 3 


10 


-1.55 


-1.48 


15.615 




2, 3, 4, 5 


15 


-1.24 


-1.47 


15.597 




2, 3, 5 


18 


-1.06 


-1.30 


15.673 




2, 3, 5 


20 


-1.51 


- 


15.611 


Bl 


2, 3 


36 


-1.43 


-1.44 


15.618 




2, 3 


39 


-1.09 




15.674 


Bl 


2, 3 


46 


-1.35 


11 


15.668 




1, 2, 3, 4 


51 


-1.34 


-1.33 


15.621 




2, 3, 5 


52 


-1.22 


- 


15.705 


Bl 


2, 3, 5 


59 


-1.06 


-1.23 


15.646 




2, 3, 5 


60 


-1.28 


-1.28 


15.520 




2, 3 


61 


-1.33 


- 


15.641 


Bl 


2, 3, 5 


62 


-1.28 


-1.14 


15.620 




2, 3, 5 


63 


-0.83 


- 


15.661 


Bl 


2, 3, 5 


66 


-1.05 


- 


15.607 


Bl 


1, 2, 3, 4, 5 


77 


-0.97 


-1.12 


15.725 




1, 2, 3, 4 


78 


-0.93 


- 


15.554 


Bl 


1, 2, 3, 4 


81 


-1.41 


-1.43 


15.656 




2, 3 


82 


-1.59 


-1.57 


15.601 




3 


84 


-1.03 


-1.19 


15.630 




1, 2, 3, 4, 5 


90 


-1.05 


-1.36 


15.639 




2, 3,5 


93 


-1.60 


-1.37 


15.640 




2,3 


94 


-1.13 


-1.37 


15.677 




2,3 


113 


-1.48 






** 




114 


-0.92 


-1.27 


15.686 




3 


115 


-1.30 






** 


3 


116 


-1.42 


-1.26 


15.685 




2,3 


202 


-1.19 




15.524 




2,3 



[Fe/H] spccl = [Fe/H]p e n from SPS. 

Bl: no [Fe/H] pll( j( is calculated for Blazhko stars with large light 
curve modulation; 

(Bl): Blazhko stars with slight light curve changes are included 
but photometric [Fe/H] is less certain; 

* : spectroscopic [Fe/H] seems to be wrong, see the text for de- 
tails; 

** : photometric data have larger scatter; 
*** : classification is uncertain; 
- : no or poor quality light curve is available. 
1) iBakos & Jurcsikl feOOd): 2) iBenkoetal 
ICorwin & Carnevl 1200 ll) : 41 ICarretta et alJ 11991 



2003) ; 3) 
5) Kaluznv 
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Fig. 2. Spectroscopic vs. photometric [Fe/H] of 18 RRab stars 
in M3. Dotted lines are the direct and inverse linear regression 
fits to the data. The correlation between the spectroscopic and 
photometric metallicities seems to be significant. As the two 
different measures of the metallicities are basically indepen- 
dent, their correlation indicates that a metallicity range of the 
M3 variables is indeed detected. 
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Fig. 3. Histogram of correlation coefficients (r) of 10000 
Monte-Carlo simulations of 18 [Fe/H] spect and [Fe/H] p h t pairs, 
a) No metallicity spread is assumed, the standard deviations of 
the data are 0.22 and 0.13 dex, respectively, b) The data range 
over 0.34 dex in metallicity with 0.18 and 0.11 dex standard 
deviations. Arrows indicate the regression coefficient between 
the measured spectroscopic and photometric metallicities. 



Pearson correlation coefficient between the two data sets is 
0.54, with 0.021 probability value. The small number of data 
makes it difficult to estimate the statistical significance of this 
result. Therefore, Monte-Carlo simulation has been performed 
in order to decide whether the correlation between the metal- 
licities is real, or it is only an artifact due to the large errors of 
the data. 

The Monte-Carlo simulation was realized as follows. Two 
artificial data sets were generated. First, the 0.22 and 0.13 dex 
standard deviations of the [Fe/H] spec t and [Fe/H]ph ot data were 



supposed to arise merely from the scatters of the observations 
around a single metallicity value. Secondly, it was assumed that 
a 0.34 dex [Fe/H] range was real, and the standard deviations of 
the spectroscopic and photometric data were smaller, 0.18 and 
0.11 dex, respectively (these values are the residual scatters of 
the regressions, see Fig.|5J. In both cases 10000 [Fe/H] spect and 
[Fe/H] p h ot data sets with 18 elements were generated. The his- 
tograms of the correlation coefficients (r) of the simulated data 
sets are shown in Fig. [3] While the observed r = 0.54 correla- 
tion occurs marginally in case a), it has significant occurrence 
frequency in case b). We can therefore draw the conclusion that 
in our sample the r = 0.54 correlation is indeed significant. 
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Fig. 4. Absolute magnitudes versus spectroscopic 

and photometric metal li cities of M3 RRab stars. 

M v = (V)int - 15.12 mag jHarrisI Il996l). Lines ind icate 
one of the theoretical ZAHB relations dCaloi et alJfl997l) . For 
the [Fe/H] p hot values at the different metallicities the faintest 
stars (which are likely on the ZAHB) indicate a similar trend 
of My with metallicity as expected, if the metallicity range is 
real. The brightest star V60 has the longest period (0.77 d) and 
possibly is already in an evolved stage of the HB evolution. 
The second brightest among the spectroscopic sample is V202 
which RR Lyrae classification is uncertain. The less convinc- 
ing evidence of an My - [Fe/H] relation for the spectroscopic 
values is possibly due to the larger errors in the [Fe/H] spect 
values. 

As the two metallicity determinations are completely inde- 
pendent, the fact that a correlation between these quantities is 
found points to the reality of the observed metallicity spread. 

3. M v - [Fe/H] relation 

The metallicity dependence of the absolute magnitudes of hor- 
izontal branch (HB) stars is well known, and is proven both 
from theoretical and observational points of view. Although the 
uniqueness, the slope and the zero point of this relation are still 
controversial, the reality of the detected metallicity spread can 
be justified by checking the connection between the measured 
absolute magnitudes and metallicities of the variables. 

A strict connection between the luminosities and metallic- 
ities of HB stars is only valid for zero age HB (ZAHB) stars, 
as the luminosity of stars evolving off the ZAHB sequence in- 
creases, and this results a true internal scatter of any measured 
My - [Fe/H] plot. The lower boundary of these plots have to 
follow the My - [Fe/H] relation valid for ZAHB stars within 
the limits of the observational uncertainties. 

It is stressed that the brightness of the stars has no influ- 
ence on either of the [Fe/H] determinations, which rules out 
any a priori connection between the values. If the metallicity 
spread of the variables is real, than the faintest magnitudes at 
the different metallicities have to follow the global trend valid 
for ZAHB stars. 



In Fig. 0] the absolute magnitudes of the RRab sample 
are plotted versus their spectroscopic and photometric metal- 
licities. For the [Fe/H] p hot values, the faintest stars at differ- 
ent metallicities follow the trend predicted by ZAHB models 
closely. Selecting ZAHB stars from the photometric sample by 
taking into account only the two faintest stars of each 0. 1 dex 
[Fe/H] bins a correlation with 0.84 correlation coefficient and a 
0.20 slope can be yield. This slope agrees extremely well with 
the model predictions which typically yield slopes within the 
0.18-0.26 range. 

The less convincing result for the spectroscopic metallic- 
ities is possibly due to the larger errors in the spectroscopic 
data. 

This result provides further evidence that a metallicity 
spread of about 0.3 - 0.4 dex of the M3 RRab stars has been in 
fact detected. 



4. Conclusion 

The statistically significant correlation between the spectro- 
scopic and photometric metallicities of the RRab variables in 
M3, which have ~ 0.6, and 0.4 dex ranges, respectively, leads 
to the conclusion that a true spread in the metallicities of the 
variables exists. This is in contrast with the strictly homoge- 
neous iron content of the cluster red giants, though, it is also 
worth noting that in an early study Cal overabundance of about 
0.6 dex of a cluster giant (vZ1397) was found JCohenl fT978i. 
This result has been neither confirmed nor rejected still today. 
Although some systematic difference between the metallicity 
scales of the red giants and the variables may exist, the com- 
parison of the observed metallicities indicates that the variables 
show a metallicity range spanning to about 0.3 - 0.4 dex more 
metal rich values. 

Spectroscopic observations and chemical composition de- 
terminations of cluster variables, apart from some AS measure- 
ments showing large scatter most probably due to observational 
errors, exist only of two clusters, to Cen and M3. In u> Cen as 
the iron abundances of the red giants have larger than 1 .0 dex 
range, the result that the RR Lyrae metallicities are not homo- 
geneous either is not surprising. It is important to repeat, how- 
ever, that no spectroscopic observations of any other globular 
cluster v ariables exist. The photo metric metallicities as com- 
piled in Kov acs & Walker! J200ll) show unaccountable large 
scatter in many clusters where no metallicity spread of red gi- 
ants is observed. This warns that a similar result, an observable 
metallicity spread of the variables might possibly occur in other 
clusters as well. 



Table 3. Mean meatallicity values of M3 variables within 
smaller and larger radial distances 





<[Fe/H] spcct > 


a 


N 


<[Fe/H] phot > 


(T 


N 


R < 2'.5 


-1.08 


0.23 


3 


-1.11 


0.01 


2 


R < 3:5 


-1.13 


0.18 


8 


-1.24 


0.14 


5 


R > 3:5 


-1.28 


0.22 


20 


-1.36 


0.13 


14 


all 


-1.24 


0.22 


28 


-1.33 


0.14 


19 
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ICorwin & Carnevl J200ll) have found that there are 12 RR 
Lyrae stars in the inner 81" region lying 0.1 - 0.2 mag below 
the ZAHB. They excluded the possibility that these stars were 
too faint because of photometric inaccuracies, however, they 
could not give any plausible reason for their low luminosities. 
As shown in Fig. @] metallicity effect may partly explain the 
faintness of these stars. The variables discussed in the present 
paper although do not show clear evidence of correlation be- 
tween their metallicities and radial distances from the cluster 
center, the innermost stars have systematically larger metallic- 
ities than the outer variables as is shown in Table[3] 

Radial differences of the HB type and of the frequency 
of blue straggler stars (BSS) have been also detected in M3, 
namely the horizontal branch of the inne r region is bluer and 
its BSS populatio n is more numerous ( Catel aiTet alJ 2001; 
iFerraro et alJll997l) . Though these phenomena still lack expla- 
nation, all of them indicate that there is a dynamical segrega- 
tion within the clu ster. Fitting the luminosity function of M3, 
iRood et ail fl999) have also found that the lower mass stars 
(M < 0.8M©; stars below the turnoff point) do not have a uni- 
form radial distribution as well. If the stellar composition, and 
central density of the stars nearer to the cluster's center are in 
some degree different, this might explain our finding, suppos- 
ing that all the spectroscopically measured cluster red giants 
with homogeneous metallicities belong to the outskirt of the 
cluster. 

If the iron content is indeed not homogeneous in HB stars, 
its explanation could probably be found in the early star for- 
mation and chemical evolution history of the system. Recent 
investigations indicate that self-pollution could have been effi- 
cient during the first some hundred million years of the clus- 
ter lifetime, just before its first crossing of the galactic plane, 
which removed all the intracluster gas from the cluster. As a 
manifestation of self-pollution, both the then existing stars ac- 
creted processed material, and the composition of the intra clus- 
ter gas which formed to sta rs later, was altered JThoul et alJ 
l2002t ID' Antona et alJEool . Models show that in M3, 11% 
of the mass of a 1 Mq star may originate from the mass loss 
of the first generation, large mass, evolved stars, and not from 
the uncon taminated protocluster gas according to the accretion 
scenario JThoul et alJ2002h . lShustov & WiebeN2000l) modeled 
supernova explosions in globular clusters, and concluded that 
self-pollution via SN ejetion could also happen in globular 
clusters, that might be a clue to explain any detected metal- 
licity spread. They also argue that SN self -pollution results in 
smaller scatter of the observed iron abundances than in any of 
the lighter elements. This helps to interpret why [Fe/H] values 
seem to be homogeneous in most of the globular clusters. 

The question why [Fe/H] dispersion only of RR Lyrae stars 
is detected, while red giants seem to be uniform in their iron 
content is difficult to answer on the bases of the presently avail- 
able observational data. Because of the small sample size and 
the large uncertainties especially in the spectroscopic data of 
RR Lyrae stars, good quality detailed spectroscopic measure- 
ments of a much larger RR Lyrae sample is needed to confirm 
and to find the correct interpretation of the suggested metallic- 
ity dispersion of RR Lyrae stars in M3. Similar studies of other 
globular clusters are also encouraged. 
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